Is there a Phase Transition to the Flux Lattice State? 
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The sharp drops in the resistance and magnetization which are usually attributed to a phase 
transition from the vortex liquid state to a crystal state are explained instead as a crossover between 
three and two dimensional behavior, which occurs when the phase coherence length in the liquid 
becomes comparable to the sample thickness. Estimates of the width of the crossover region and 
the phase coherence length scales are in agreement with experiment. 
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For some years now it has been widely believed that the 
vortex liquid state of a type-II superconductor changes 
on cooling to a flux lattice at a first-order phase transi- 
tion [pi . Sharp drops have been seen in the magnetiza- 
tion |]-||] , entropy Q and resistance |(| and it is usually 
assumed that such drops are due to a freezing transi- 
tion. I shall describe here a simple alternative expla- 
nation of these drops, based on the idea that they are 
caused by a crossover from three-dimensional behavior 
to two-dimensional behavior when the phase correlation 
length along the field direction in the vortex liquid, l\\ , be- 
comes comparable to the sample dimensions |7j . Because 
l\\ grows very rapidly as the temperature is lowered the 
crossover region appears narrow enough to masquerade 
as a first order phase transition. An estimate of its width 
is given and found to be comparable with the width of 
the drop in the magnetization || . 

Central to my picture is the idea that at all non-zero 
temperatures the system is in the vortex liquid state and 
that it is only finite size effects which produce the sem- 
blance of a phase transition. In the thermodynamic limit 
of an infinitely big system there would be no phase tran- 
sition, no sharp drops and only a vortex liquid phase. 
At zero temperature in the ground state there would of 
course be the usual Abrikosov lattice. The length scale 
over which there is short-range crystalline order in the 
vortex liquid, l±, is expected to increase without limit as 
the temperature is lowered towards absolute zero. 

Essentially the picture being advocated here is that 
of Refs. |,|, where it was shown that in both two and 
three dimensions phase coherence in the vortex crystal 
state was destroyed by the phase fluctuations associated 
with the shear modes of the crystal and it was suggested 
that as a consequence crystalline order might not exist 
in these dimensions. Since it is not obvious why the loss 
of phase coherence need necessarily imply the absence 
of crystalline order this suggestion has (to say the least) 
not been widely accepted. However, our Monte Carlo 
simulations in two-dimensions done within the 

lowest Landau level approximation (LLL), found no fi- 
nite temperature phase transition and demonstrated that 
the length scale over which phase coherence is lost is l±, 



which diverged in the zero-temperature limit. 

Strong support for the idea that the apparent first- 
order phase transition is just a crossover from three 
to two-dimensional behavior comes from multi-terminal 
transport (flux transformer) measurements on untwinned 
single crystals of YBCO |j and on thin BSCCO crys- 
tals juj. In these experiments current is injected and 
removed at two electrodes attached to the top face of 
the crystal and voltage differences measured between 
two points lying between the electrodes on the top face, 
Vtop, and between two points directly below them on 
the bottom face of the crystal, Vbot- The magnetic field 
is along the c-axis of the crystal and the thicknesses of the 
YBCO crystals ranged up to 35 [im while the BSCCO 
crystals were around 1 /im thick. For crystals of these 
thicknesses both voltages became very nearly equal to 
each other as the temperature was reduced after which 
both fell sharply to zero. The vanishing of the voltage is 
usually taken as marking the onset of the freezing transi- 
tion into the crystalline state. When Vtop and Vbot are 
equal it means the flux lines are moving without cutting 
and reconnecting so there must be phase coherence across 
the crystal. As Vtop an d Vbot are becoming equal as 
the "freezing" transition is approached, then the length 
scale over which there is phase coherence in the direction 
of the magnetic field, l\\, must be growing as the temper- 
ature of the vortex liquid is reduced, becoming compara- 
ble to the crystal thickness right at the apparent phase 
transition. Furthermore, measurements of the non-linear 
current-voltage characteristics on the same crystals give 
information on the length scale l± via the scaling rela- 
tion aE = Jf(Jl±l\\&o/kBT), where E is the electric 
field perpendicular to the field, J is the current density 
in the same direction, a is the conductivity and $o is 
the flux quantum [lq] . Preliminary studies indicate that 
it too is a length of micron magnitude p6| . In the con- 
ventional picture of a first- order phase transition it is 
very hard to understand why long-length scales should 
be present in the vortex liquid — /|| is at least as large 
as 35 \xm in YBCO at the freezing transition — and why 
there should be a coincidence of the freezing transition 
with the temperature at which this length scale becomes 
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comparable to the dimensions of the sample. My hy- 
pothesis that the "freezing" phenomenon is a crossover 
induced effect seems much more natural. 

We start from the Ginzburg-Landau free energy for a 
superconductor in a uniform magnetic field B = V x A; 

F[¥] = J d 3 r(a(T)|Mf + ||*| 4 

(-ihd M - 2eA^\ 2 
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Here a(T) is a temperature-dependent variable e.g 
a(T) — — a'(l — T/Tco), T c q is the mean- field transi- 
tion temperature, j3 is the coupling constant and m M is 
an effective mass. The masses in the ab plane will be 
taken equal and denoted by m a f, and the mass in the c 
direction is m c . To start with we shall neglect the fluc- 
tuations in the vector potential A and restrict the order 
parameter to the space spanned by the lowest Lan- 
dau levels (LLL). Then aji = ol(T) + eKB/m a b = de- 
fines the mean- field line H C 2. The temperature is conve- 
niently and conventionally JTt) represented by the dimen- 
sionless parameter ar = a B m r ) 1 ■ 

Low temperature is represented by ax — > — oo, high 
temperature by ax — > oo. More conveniently ar — 
(2/Gi) 1 / 3 (ht)- 2 / 3 (l - t-h) where Gi is the Ginzburg 
number, (about 0.02 for YBCO fjj), t is the reduced tem- 
perature T/T c q and h = B/B C 2(0) where B C 2(0) is the 
straight line extrapolant of the B c i to zero temperature 
(about 168 T in YBCO §). 

The length scale 2ii is a measure of the phase 
correlation length along the field direction (and is 
technically the length scale on which the propagator 
{9*{x, y, z)ty(x, y, z')) decreases as a function of \z — z'\. 
The other important length scale in the problem, 
which measures the degree of short-range crystalline or- 
der in the vortex liquid can be obtained from the density- 
density correlation function (^(x, y, z)\ 2 \^(x', y', z)\ 2 } c . 
In addition l±_ seems to be equal to an (appropriately de- 
fined) phase coherence length scale in the direction per- 
pendicular to the field [Q . It is our contention that these 
lengths reach infinity as the temperature goes to zero and 
that as a consequence it is useful to use the techniques 
of "zero-temperature scaling" jL8| . 

At zero-temperature the ground state order parame- 
ter $!o(x,y) which minimises F[^f] corresponds to the 
Abrikosov flux lattice of straight vortex lines arranged in 
the form of a triangular lattice. Expanding about the 
minimum by setting V I / = V&o + S^f[x, y, z] one finds ^|Jl9| 
that the low energy states above the minimum can be 
described by an effective Hamiltonian 



H = \J d 3 v(^p s {d6/dz) 2 + c 66 / 4 (V 2 



(2) 



Here 9 is the phase change associated with S^f, p s is the 
superfluid density (p s = p s h 2 \aH\/rn c P in the LLL mean- 



field limit, with p~ s = I/Pa and Pa = 1.1596...) and 
cqq is the shear modulus equal in LLL mean-field ap- 
proximation to 2cQQa 2 H /f3 with cqq = 0.0885... . The 
magnetic length scale I = y/ft/2eB. The relation be- 
tween the phase change 9 and flux line displacements 
is u x = —l 2 d9/dy and u y = l 2 39/dx ||. I believe 
that vortex crystals behave differently from other crys- 
tals because only for them can displacements be written 
as phase derivatives. 

It is convenient to intr oduce the dimcnsionless vari- 
ables Z where z = £ C Z; £ c = (ft 2 /2m c \aH l) 1 ^ 2 and 
R = (X,Y) where x = IX, y = IY . Then the effective 
Hamiltonian in d dimensions becomes 



k B T 2T 



d d - 2 T,\p s (d9ldZ^ 



(3) 



The subscript i takes values 3,4, ... , g? and an implied 
summation convention is being used. In three dimen- 
sions the dimensionless effective temperature T equals 
4ir/\a T \ 3/2 . 

For the case of a thin film of thickness L z where the 
variation of ^ with respect to z can be ignored, the 
Hamiltonian reduces to 



H 



1 



k B T 2T 



d 2 Rc 66 {V 2 ± t 



(4) 



where in two dimensions the effective temperature T 
equals 7r/a2T> with olit — olh{^L z / peBksT) 1 ^ 2 . 

Hamiltonians such as Eq. (^) have been studied in con- 
nection with smectic liquid crystals p0| where a renor- 
malisation group treatment has proved very useful. To 
this end rescale lengths so that x' = bx,y' = by but 
z' = b 2 z when in order to keep the free energy invari- 
ant the renormalised temperature T' = b 6 ~ 2d T. Setting 
b = e p one obtains the following flow equation for the 
effective temperature as a function of the length scale: 



dT 
dp 



(6 - 2d)f. 



(5) 



For dimensions d < 3 the flow takes one away from low 
temperatures to the high-temperature sink. Above three 
dimensions T = is a stable fixed point indicating that 
3 is the lower critical dimension for behavior controlled 
by this fixed point. 

Consider first the case of thin films (d = 2). Start- 
ing at a low effective temperature T where the correla- 
tion length is l±, scale up until the effective tempera- 
ture T 1 is a, of order unity, and where the correlation 
length is I. Then b = l ± /l and a = f' = {l±/l) 2 f so 
l± ~ lil/f) 1 / 2 ~ l\a 2T \. This result for l x is consis- 
tent with our Monte Carlo simulations in two dimen- 
sions 10- 1^] and has also been derived within the so- 
called parquet graph approximation There have 
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been a number of other Monte Carlo simulations in two- 
dimensions employing (the experimentally unrealisable) 
quasi-periodic boundary conditions and the LLL approx- 
imation in which a first order transition to a crystalline 
state seems to take place at a finite temperature Q . As 
no first-order transition has been reported for supercon- 
ducting films, I shall take the apparent phase transition 
to be an artifact of the simulations. For a review of this 
controversy, see Ref. Q. 

In thin superconducting films ac resistivities |23j | and 
flux liquid viscosities [^i) have been apparently success- 
fully interpreted in terms of the continuous phase transi- 
tions predicted in conventional two-dimensional melting 
theory pjj], which would seem to be at odds with the 
idea that the vortex liquid state should be the only state 
present at all finite temperatures. These experiments es- 
sentially determine the relaxational timescale of the sys- 
tem which grows rapidly at low temperatures. According 
to our Monte Carlo simulations 1 1 2(| it increases with 
an Arrhenius form in the two-dimensional effective tem- 
perature, i.e. as exp(C/T), with C a constant of order 
one. Such an expression provides an adequate fit to the 
experimental data. Invoking a finite temperature phase 
transition thus seems unnecessary. 

Returning now to bulk systems one sees that the term 
on the right-hand side of Eq. (^J) is zero when d = 3. 
To describe three dimensions it is necessary to go be- 
yond the simple scaling flow equation. At low temper- 
atures one expects it to become instead J2?j], dT/dp = 
T 2 /(2irA) where A is a constant. On integrating up 
to the lengthscales at which T' is of order unity, where 
b = IJl = (l\\/tc) 1/2 , one has l\\ ~ £ c exp(A|a T | 3 / 2 ) 
and l±_ ~ I exp(^4|a-r| 3 / 2 /2). Such exponentially rapid 
growth of correlation lengths is typical of a system at its 
lower critical dimension. It is perhaps noteworthy that 
phase coherence length scales of this form were previously 
obtained via estimates of the crossing energy of flux lines 

When the sample thickness L z is greater than lu then 
the system behaves as a three-dimensional bulk system, 
but as the temperature is lowered hi grows and eventu- 
ally becomes greater than L z . Then the vortices behave 
as in a two-dimensional thin-film system with phase co- 
herence across the sample. The crossover region is where 
L z ~ /||. If one knew A one could get a good estimate 
of otfp, the value of olt at the crossover, for a sample of 
known thickness . At this time we shall take a* T from ex- 
periment p9fl . It is around —9 for YBCO crystals of (say) 
40 pm dimensions |30| . It follows from the definition of 
«t that the crossover line obeys the LLL scaling relation 
B ~ (T c0 — T) 3 / 2 for fixed L 2 , except for temperatures 
so close to T c q that it is no longer possible to regard £ c 
as a constant. 

Because of the exponentially rapid growth of lu, the 
width of the crossover region as a function of field and 



temperature is quite narrow. Its width may be estimated 
by finding the small change Sax in a T to make l\\ grow 
to 2L Z . Then 2>Sa T /2a* T ~ ln(2)/ ln(L 2 /£ c ) - SB/B, 
where SB is the change in B required to double lu. For 
the 0.2 mm thick crystal of Welp et al. M, we predict that 
SB/B is about 0.054 at a field of 4.2 T while the observed 
width looks to me to be around 0.2 T, making the ratio 
SB/B about 0.05. Thus the crossover width is similar to 
the rounding of the jump in the magnetization. If one be- 
lieved in a first-order transition one would have to claim 
that the rounding of the jump was somehow caused by 
a spatial inhomogeneity of oxygen concentration or some 
non-uniformity of the applied field. Another feature of 
our crossover argument is that as the sample thickness is 
increased then the crossover should occur at lower tem- 
peratures, but the decrease of the "transition" temper- 
ature is only as the logarithm of the sample thickness. 
The magnitude of the effect is easily found from solving 
L z ~ lu and certainly looks measurable provided that a 
range of crystals can be made of various thicknesses but 
otherwise identical. 

When the crossover takes place the effective tem- 
perature of the resulting two-dimensional (thin film) 
state is very low. Using the definition of U2T 
to re-express it in terms of h and t as olit = 
0.5(Gi/2) 1/12 (/rf) 1/6 a T (£ z /£o) 1/2 , where £ is the zero- 
temperature c-axis correlation length, (h /2m c a') 1 ' 2 
(which is about 2 A) and then inserting numbers ap- 
propriate to YBCO at a field of 4 T and sample thick- 
ness 0.2 mm one finds a2T ~ —1600. Notice that in 
the crossover region and below, l± is more than 10 3 I. 
It would then only require a tiny amount of coupling 
between the vortex system and the underlying real crys- 
tal to produce a system whose structure factor displayed 
Bragg spots rather than the rings expected for a liquid 
state pjj . Thus the vortex system would appear crys- 
talline to most experimental probes. 

Because the two-dimensional state below the crossover 
is at such a low effective temperature, its entropy can 
be obtained by differentiating the mean-field free-energy 
per unit volume, —a 2 H /2/3(3A, with respect to the tem- 
perature. However, the entropy of the three-dimensional 
state when olt = a* T ~ — 9 is a few percent above the 
mean- field value ]32j| . Therefore the entropy must de- 
crease rapidly on cooling through the crossover region, 
giving rise to an apparent jump. Per vortex per double 
layer the entropy change AS can be written as 



AS ~ fc s ( S /e )(2/«) 1 / 6 t 2 / 3 / l - 1 / 3 5 (a T ), 



(G) 



where h,t and ax are at their crossover values, s is 
the distance between the double layers and the function 
g(ax) has at present to be inferred from the numerical 
simulations Q of specific heats. I estimate that AS is 
around ks at 4 T. Notice that Eq. (||) predicts that AS 
increases as the field is reduced in accord with experiment 
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The argument of this paper is that the properties of 
the vortex liquid are controlled by a particular zero- 
temperature fixed point. Because in dimensions two and 
three there are diverging length scales associated with 
this fixed point, the concept of renormalisation group 
relevance and irrelevance can be used to determine the 
importance of effects ignored in the calculation. Thus 
higher Landau levels are clearly irrelevant in the low- 
temperature limit as they are associated with more mas- 
sive fields. Similarly screening effects — the fluctuations 
of the vector potential A — can be shown to be irrelevant 
at this zero-temperature fixed point by an extension of 
the argument used in Ref. |33f| , where such fluctuations 
were shown to be also irrelevant at the critical fixed point. 
Thus our treatment of the length scales l\\ and l± should 
become exact in the low-temperature limit. However, if 
there are terms present which couple the vortices to the 
underlying crystal (as in numerical simulations done on 
a grid) then a spectrum of low-energy excitations may 
arise different from that in Eq. (||) which could result in 
a finite temperature phase transition. 

One important effect so far neglected is disorder. Af- 
ter electron irradiation, which induces point defects, the 
"drops" in the resistance are replaced by a smooth de- 
crease with temperature [[mJ. Disorder might be ex- 
pected to limit the growth of l± and Z|| at the Larkin 
length scales R± and Ru [Q, which would naturally pre- 
vent a sharp crossover effect from occurring. Even in 
non-irradiated samples the resistance drop is much nar- 
rower as a function of temperature than the drop in the 
magnetization. I believe that the sharper drop in the re- 
sistance is a consequence of a percolative mechanism in 
which regions of transverse dimensions l± and longitudi- 
nal size Z|i become pinned by the disorder and when these 
pinned regions percolate across the system then the resis- 
tance falls away. In thicker (90 /j,m) untwinned crystals 
of YBCO Vtop does not equal Vbot when the resistance 
drops to zero |35|] — a fact consistent with a percolative 
mechanism. 

Thus it seems possible to give an account of the main 
features of the mixed phase at a semi-quantitative level 
by employing the ideas of zero- temperature scaling. 
Should there actually be a finite temperature phase tran- 
sition then these arguments become irrelevant. However, 
it is hard to see how the long length scales revealed by the 
flux transformer experiments can be incorporated into 
the conventional melting picture. 
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